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We study the dynamics of a single photon pulse travels through a linear atomic chain coupled to 
a one-dimensional (ID) single mode photonic waveguide. We derive a time-dependent dynamical 
theory for this collective many-body system which allows us to study the real time evolution of 
the photon transport and the atomic excitations. Our analytical result is consistent with previous 
numerical calculations when there is only one atom. For an atomic chain, the collective interaction 
between the atoms mediated by the waveguide mode can significantly change the dynamics of the 
system. The reflectivity of a photon can be tuned by changing the ratio of coupling strength and the 
photon linewidth or by changing the number of atoms in the chain. The reflectivity of a single photon 
pulse with finite bandwidth can even approach 100%. The spectrum of the reflected and transmitted 
photon can also be significantly different from the single atom case. Many interesting physical 
phenomena can occur in this system such as the photonic bandgap effects, quantum entanglement 
generation, Fano-like interference, and superradiant effects. For engineering, this system may serve 
as a single photon frequency filter, single photon modulation and may find important applications 
in quantum information. 

PACS numbers: 42.50.Nn, 42.50.Ct, 32.70.Jz 


I. INTRODUCTION 

Photons are ideal carriers for quantum information. 
Hence, manipulating and routing photons can have im¬ 
portant applications in the quantum information tech¬ 
nology which has been extensively studied M- How¬ 
ever, the photons rarely interact with each other. There¬ 
fore, we have to resort to the photon-atom interaction 
to modulate and control the photons. The photon-atom 
coupling in the vacuum is usually very weak. We can, 
however, modify this coupling strength by changing the 
environment of the vacuums by Purcell effect Q. Strong 
photon-atom interaction can be achieved by confining the 
photon in reduced dimensions such as in quasi ID pho¬ 
tonic waveguide with transverse cross sections on the or¬ 
der of a wavelength square |^. There are a number 
of systems that can act as a quasi ID waveguide such 
as optical nanofibers [l^, photonic crystal with line de¬ 
fects [m, surface plasmon nanowire fl3. and supercon¬ 
ducting microwave transmission lines This kind 

of waveguide is extremely interesting because it can not 
only enhance the interaction but also guide the photon 
which is important for the information transport. 

A single photon scattered by a single atom embedded 
in a one-dimensional waveguide has been studied in the 
pioneering work by Shen and Fan [il,[i3. They em¬ 
ployed the real space description of the Dicke Hamilto- 
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nian and Bethe-ansatz to derive the stationary transport 
properties of a single photon-atom interaction. It shows 
that the photon transmission spectra can be strongly 
modified. In particular, a photon with frequency reso¬ 
nant to the two-level atom can be even completely re¬ 
flected. Since then, this method has been generalized to 
the case for multi-level atom and multi-photon interac¬ 
tions [TtI - I^ . In the stationary calculations, the pho¬ 
ton is assumed to be a plane wave with single frequency. 
However, in the real case the photon is always a pulse 
with a finite bandwdith. An extension of the station¬ 
ary theory was developed by Rephaeli et al. to derive 
the perfect inverting pulse where they project the initial 
quantum state to the stationary eigenstates of the sys¬ 
tem and then apply a time evolution operator In 

2011, Chen et al. illustrated a time-dependent theory 
for a single photon absorption by single emitter coupled 
to one-dimensional photonic waveguide [^ . This the¬ 
ory can be used to study the dynamics of the photon 
wavepacket and the atomic excitation in real time. The 
spontaneous emission from a pair of two-level atoms near 
a nanofiber has also been studied where a substantial ra¬ 
diative exchange is predicted between distant atoms [1^ . 

The collective spontaneous emission by an ensemble of 
identical atoms has been extensively studied since the 
pioneering work of Dicke in 1954 (26l - K^ . The Dicke 
super-radiance provides a valuable example for studying 
the many-body physics of photons and atoms [s^- The 
interaction between the atoms mediated by the common 
vacuum field can dramatically change the d ecay rate and 
the eigenspectrum of the atomic ensemble [35l| . The de¬ 
cay rate of the system can be superradiant or subradiant. 
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FIG. 1: (Color online) Single photon transport through a 
linear atomic chain coupled to a one-dimensional photonic 
waveguide. The atoms can be coupled by the waveguide vac¬ 
uum modes. 


The Dicke symmetric state of maximum cooperation has 
decay rate N times the single atom decay rate [1^. The 
emission spectrum can also be significantly altered. The 
stationary calculation of a single photon scattered by 
multiple emitters coupled to one-dimensional photonic 
waveguide has also been studied where it is shown how 
the transmission spectrum can be strongly modified by 
the quantum interference effects [3^1^ . However, in this 
analysis the collective dipole-dipole interactions between 
the atoms are not included. 

A general master equation for describing the interac¬ 
tion between a quantum system and a continuous-mode 
Fock state has been studied [S^. In this paper, we use 
a different approach to derive a time-dependent theory 
of a single photon scattered by an atomic chain coupled 
to a ID waveguide where the collective many-body in¬ 
teraction can play a crucial role. We can calculate how 
the atomic excitations evolves in real time and how the 
single photon pulse is reflected and transmitted. We also 
show how the reflectivity and the spectrum are modi¬ 
fied by this collective many-body system. Many inter¬ 
esting physical phenomena can occur in this system such 
as the photonic bandgap effects, quantum entanglement 
generation, Fano-like interference, and superradiant ef¬ 
fects. This system may serve as a single photon frequency 
filter, single photon modulation and may find important 
applications in quantum information. 

This paper is organized as follows. In section II, we 
derive the time-dependent theory for a single photon 
transport through an atomic chain coupled to a one¬ 
dimensional photonic waveguide. In section III, we dis¬ 
cuss the one-atom case where we derive an analytical 
solution which has not yet been obtained before and our 
result is consistent with the previous numerical results. 
In section III, we numerically study how a single photon 
can be scattered by a two-atom system and show some 
new features compared to one-atom case. In section IV, 
we study the N-atom case and show how the atomic sep¬ 
aration and atom numbers affect the reflection and trans¬ 
mission of a single photon. In section V, we study the 
effects of spontaneous decay to the free space. Finally, 
we summarize the results. 


II. THEORY 

We consider a linear atomic chain coupled to a one¬ 
dimensional single-mode waveguide (Fig. 1). We assume 


that the atoms are equally spaced with separation a and 
the number of atoms is Na- A linearly polarized single 
photon pulse with spectral width Aw propagates through 
and interacts with the atomic chain. The Hamiltonian of 
the system in the rotating wave approximation is given 
by Kii 

Na 

H = h{oja-i]r)'^S^S~ + ny^utkajak 

j^l k 

Na 

+^E (1) 

j=l k 

where the first term is the atomic energy including the 
spontaneous decay to the free space, the second term 
is the energy of the guided photon, the third term is 
the coupling between the atoms and the guided photon 
modes. Here Wa is the transition frequency of the two- 
level atoms, 7 is the rate of spontaneous decay to the non- 
guided modes, = \e)j{g\{S~ = \g)j{e\) is the raising 

(lowering) operator of the jth atom, aj,(a^) is the cre¬ 
ation (annihilation) operator of a photon with wavevec- 
tor k, gk is the coupling strength between the atoms and 
the guided photon mode. Wfc is the angular frequency of 
the photon with wavevector k and rj is the position of 
the jth atom. Assuming that Wa is far away from the 
cutoff frequency of the photonic waveguide and guided 
photon has a narrow bandwidth, the dispersion relation 
for the guided photon can be approximately linearized 
as w/c = Wa + (IA:I — ka)vg where ka is the wavevector at 
frequency Wa and Vg is the group velocity [40l| . 

For the single photon excitation, the quantum state of 
the system at any time can be expressed as 

No. 

|d/(t)) = ^a,(t)e—“*|e„0)+^/3,.(t)e-“'"‘|g,U) 

j=l k 

( 2 ) 

where |ej,0) is the state in which all the atoms except 
the jth atom are in the ground state with zero photon 
in the waveguide and the free space, \g, Ik) is the state 
in which all the atoms are in the ground state and one 
photon with wavevector k is in the waveguide. 

The dynamics of the system is given by 

(3) 

k 

Na 

Mt) = (4) 

1=1 

where 6uJk = Wfe — Wa = (|fc| — ka)vg is the detuning be¬ 
tween the atomic transition frequency and the frequency 
of the guided photon. 

Integrating Eq. (4) we obtain 

Na 

Pk{t) = fik{Q)-iy^gle-^’^^^ / , (5) 

i=i -^0 
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where /3fe(0) is the initial photon amplitude. On substi¬ 
tuting from Eq. (5) into Eq. (3) we get the dynamic of 
the atomic excitations 


where Sk = |fc| — ka with k < 0. 

To calculate the photon spectrum when t —>■ oo, we 
define 


aj{t) 


k 

pOO 

EE \g^\2e-^Hr,-rO / 

*=i ° then have 


/ OO 

( 12 ) 

-OO 

= 0 when t < 0 with j = 1, 2, • • • , Na- We 


( 6 ) 


For a long ID waveguide we can replace the summation 
over k by integration 




Xj{5k)e-^^'^'’o^d5k. 


(13) 



(7) 


Inserting Eq. (13) into Eq. ( 8 ) we obtain a set of linear 
equations which are given by 


where L is the quantization length in the propagation 
direction. 

According to the Weisskopf-Wigner approximation, on 
summing over k and integrating over t', Eq. ( 6 ) becomes 

No. j, 

= hAt) - E( 2 "*'“’'^' + 

1=1 ® 

with 

hAt) = -i\/^ (9) 

being the excitation by the input photon, E = 
2L\gkJ‘^/vg being the coupling strength between the 
atom and the guided photon and rji = \rj — ri\. The 
second term is the collective many-body coupling be¬ 
tween the atoms induced by the guided photon modes. 
It is noted that the collective interactions induced by the 
guided photon is a long-range effect and there is a time 
retarded effect for these interactions. The real parts of 
the collective coupling give the collective damping while 
the imaginary parts yield an energy shift. Given the ini¬ 
tial conditions, we can calculate the atomic excitations 
at an arbitrary time from Eq. ( 8 ). However, the analyt¬ 
ical solution for Eq. ( 8 ) is usually difficult and we need 
to resort to numerical methods such as Finite-Difference 
Time-Domain (FDTD) [4l|. 

Knowing the atomic excitation aj(t), we can also cal¬ 
culate the amplitude of each guided photon mode at an 
arbitrary time using Eq. (5). To be specific, if the input 
photon comes from the left end, the right propagating 
modes at time t are given by 

I3g{t) = / 3 ,( 0 )-*y^ ^-dk.+Sk)r, 

( 10 ) 

where Sk = k — ka with fc > 0. For the left propagating 
modes we have 





( 11 ) 


Na 


- iSkvgXjiSk) = bj{6k) - ^Xji^k) + ^ VA[{Sk)xi{Sk) 

1^1 

(14) 

where 


bASk) = -i (15) 


and the collective coupling matrix V is given by 


VNASk) = —, 


1 

Aka 


^ika ^ ^ ^ A(Na — i-)ka 


1 


^i(Na — l)ka ^i{Na—‘2)ka 


^i(Na—2)ka 


1 


(16) 

with k = ka + Sk. 

Let M{Sk) = —VnAA) + ( 7/2 — iSkvg)lNa where 
is a Na X Na unit matrix. We have 

Na 

xASk)^^M-ASk)biiSk) (17) 

1=1 

We can then calculate the photon spectra for the right 
and left propagating field at t —>■ 00 : 


^ 1=1 

(18) 


(19) 

1=1 


We can calculate the photon pulse shape after the in¬ 
teraction by the Fourier transformations as given by 



/ OO 

(ik{t)A^^^"^'’^'^dSk. ( 21 ) 

-OO 
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We can also calculate the reflectivity R and the trans- 
mittivity T of this system and they are given by 


T f°° 

\l3h{t^^)\^d6k, ( 22 ) 

r poo 

\l3B.{t^^)\^ddk. (23) 


In a pure ID photonic waveguide or a quasi-lD pho¬ 
tonic waveguide with F ^ 7 , we can neglect the sponta¬ 
neous decay ( 7 ) to the non-guided modes. In the follow¬ 
ing sections, we first neglect 7 and mainly focus on the 
collective effects induced by the guided modes. We ana¬ 
lytically and numerically study the dynamics of a single 
photon pulse interacting with one atom, two atoms and 
N equally spaced atoms. In Sec. VI, we briefly discuss 
the effects of the spontaneous decay to the free space. In 
all calculations, we assume that the atoms are initially in 
the ground state and the input single photon pulse has a 
Gaussian shape which is given by 


^5k{0) = 


(87r)i/4 


{ka + 5k- 
- 



(24) 


where fep is the wavevector corresponding to the center 
frequency of the pulse and A is the width in the k space 
with the full width at half maximum of the spectrum 
being ^2 ln2Ai;g. A single photon with Gaussian pulse 
shape can be generated by single photon electric-optic 
modulation as has been experimentally demonstrated in 
[ 4 ^ . We note that ^ \Psk{0)\‘^dSk = 1 which is the 

requirement for a single photon number. 


III. SINGLE ATOM 

In [^, Chen et al. have numerically studied the case 
of a single Gaussian photon pulse scattered by a single 
atom coupled to a ID waveguide. Here we present some 
analytical results based on the method shown in the pre¬ 
vious section. 

From Eq. ( 8 ), the atomic excitation for one-atom case 
with 7 = 0 is governed by 


ai{t) = bi{t) (25) 

where bi{t) is given by Eq. (9) with j = I. The formal 
solution of this equation is given by 

a{t) = —zS'{erf(C' -I- Vsto) — erf(C')}e“ 2 * (26) 

where C = {A — 2Bto)/2'/B with fg = 

A = —r/2 -I- i{ko — ka)vg, and B = A^Ug/4. 

S = , and erf(a;) = 

(2/y^) Jq e~* dt is the error function. This result has 
been derived in Ref. with ko = ka- Here we gener¬ 
alize this result to include the case when ko ^ ka, i.e.. 


when the center frequency of the input photon pulse is 
not equal to the transition frequency of the atom. 

The atomic excitation as a function of time when T = 
Avg and fcg = ka is shown as the black solid line in Eig. 
2(a). When t = 0, the center of the input photon pulse 
is 10/A away from the atom. The atom is first excited 
and then deexcited when the photon pulse propagates 
through the atomic chain. The photon pulse incident 
from one direction can be separated as odd and even 
modes [ 2 ^ . Since only even mode can couple to the atom, 
the maximum excitation by a photon incident from one 
direction is 50%. The maximum atomic excitation for 
the Gaussian input in our example is about 40% which is 
close to the 50% limit. The 50% limit can be reached by 
the inversion pulse incident from one side as discussed in 
Ref. [ 2 ^ (see the red dotted line in Fig. 2(a)). 

For one atom, Eq. (14) becomes 

- iSkvgXi{Sk) = bi{Sk)-^Xi{dk) (27) 
whose solution is 


Xi{Sk) = -i 



fdskjO) 

^ — iSkvg 


(28) 


The photon spectra after the interaction can then be cal¬ 
culated from Eqs. (18) and (19), and given by 


00 ) 

00 ) 


{2iSkvg/T)l3sk{0) , , 

1 - 2tSkvg/r ’ ^ ’ 

2i{ka+Sk)ri I^Sk ( 0 ) _ 

l-2i6kvg/T' ^ ’ 


with (3sk{0) given by Eq. (24). The spectra of the pho¬ 
ton before and after the interaction when /cg = and 
F = Avg are shown in Fig. 2(b) where we see that the 
reflected spectrum is similar to the initial spectrum but 
the transmitted spectrum has two peaks. However, the 
reflected spectrum is narrower than the incoming spec¬ 
trum and the reflected bandwidth can be controlled by 
changing the coupling strength. The smaller the cou¬ 
pling strength becomes, the narrower the reflected band¬ 
width is. In this case, there is about 66 % possibility 
that the photon is reflected. When Sk = 0, we have 
oo)P = l^5fe(0)P and |/3||(t oo)p = 0 which 
means that the resonant frequency component is com¬ 
pletely reflected and it is independent of the input pulse 
shape. This is consistent with the stationary calculations 
in Ref. [l^. The positions of the peaks 6k± in the trans¬ 
mitted spectrum are given by 


Sk± _ ^ T] 
A ~ 2V2 




(31) 


where 77 = F /Avg is the the ratio between the coupling 
strength and the input photon bandwidth. We can see 
that the peaks depends only on this ratio. When 77 —>■ 0, 
Sk± —>■ 0. The splitting of these two peaks increases 
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FIG. 2: (Color online) (a) The atomic excitation as a function time where black solid curve is the result for Gaussian input 
while the red dotted curve is the result for the inversion pulse given by Ref. [^. When t — Q the center of the input pulse is 
10/A away from the atom, (b) The spectrum of the incoming, reflected and transmitted photon when F = Aug and ko = ka. 
(c) The pulse shapes of the incoming (at t = O/F), reflected (at t = 20/r) and transmitted (at t = 20/r) photon when F = Avg 
and ko = ka- The atom is at the position x = 0. (d) The reflectivity, transmittivity and relative full width at half maximum 
of the reflected pulse as a function of coupling strength. The unit of the reflection bandwidth is \/2 In 2Avg. (e) The spectrum 
of the incoming, reflected and transmitted photon when F = Avg and ko — ka = A/2, (f) The reflectivity and transmittivity 
as a function of detuning ko — ka when F = Avg. Here, I is incoming; R is reflection; T is transmission. 


when we increase 77 and it approaches when 77 —^ 00 . 
When f] = 1, 6k± = which are the peak positions 
of the transmitted spectrum shown in Fig. 2(b). 

The photon pulse shape before and after the interac¬ 
tion can be calculated by Eqs. (19) and (20) and they 
are shown in Fig. 2(c) where the black solid curve is the 
incoming photon pulse while the red dotted curve is the 
reflected photon pulse and the blue dashed curve is the 
transmitted photon pulse. The reflected pulse shape is 
similar to the incoming pulse shape while the transmitted 
pulse shape is quite different where there are two peaks 
due to the interference between the incoming photon am¬ 
plitude and the reemitted photon amplitude. 

The reflectivity and transmittivity of this system, ob¬ 
tained from Eqs. (21) and (22), are given by 


R = 



e-^y" 

1 -b V/??^ 


dy, 


T = 1-R. 


(32) 

(33) 


The reflectivity and transmittivity of a single atom sys¬ 
tem also depends only on the ratio F / Avg. When 77 —?> 0, 
i? —>■ 0 and T —>■ 1, i.e, the photon can completely trans¬ 
mit when there is no coupling. When 77 —>■ 00 , i? —>• 1 
and T —> 0, i.e, the photon is completely reflected when 
the coupling strength is very large comparing with the 
photon bandwidth. The reflectivity and transmittivity, 
as a function of the ratio T/Avg, are shown in Fig. 
2 (d) where we see that the reflectivity increases when 
F/Awg increases and it can approach 100 % reflectivity 


when F /Avg ^ 1. The reflection bandwidth can also 
be controlled by changing the ratio T/Avg. The smaller 
the F /Avg is, the narrower the reflection bandwidth is 
(see the green line with triangle symbol in Fig. 2(d)). 
This may be used to generate a single photon with ultra¬ 
narrow linewidth. 

We also study the case when the center frequency of the 
incoming photon is detuned from the resonant frequency, 
i.e., fco ka- For example, the result when kg — ka = A/2 
is shown in Fig. 2(e) where we can see that the spectrum 
is quit different. Both the reflected and transmitted spec¬ 
tra are asymmetric. However, the resonant frequency 
component here is also completely reflected. The reflec¬ 
tivity and transmittvity, as a function of center frequency 
detuning fco ~ when F = Avg, are shown in Fig. 2(f) 
where we see that the reflectivity decreases as the detun¬ 
ing increases and the highest reflectivity here is about 
66 %. The full width at half maximum of the reflection 
curve is about l.SAug. 


IV. TWO ATOM 


In this section, we study how a single photon pulse is 
scattered by a two-atom system coupled to a ID waveg¬ 
uide. This is the simplest collective many-body system. 

From Eq. ( 8 ), the atomic excitation dynamics for two- 





















6 


atom system are given by 

d{t) = b{t) + V2{Sk = 0)d{t — —), (34) 

Vg 


where d = (ai,a 2 )^, b = and the effective 

collective coupling matrix is 


V 2 iSk = 0) = -- 


1 

^ikaCL Y 


(35) 


with a being the atomic separation. The analytical solu¬ 
tions for the coupled equations (33) are difficult to obtain. 
We can apply numerical method to calculate the atomic 
excitation in arbitrary time. 

To calculate the spectrum, we follow the procedures 
shown in Sec. II. For two-atom system, Eq. (14) becomes 


— i5kvgx{5k) = b{Sk) + V2{Sk)x{5k) (36) 


where x = (xii X 2 )^, b = ( 6 i, 62 )^, and V 2 { 6 k) is given by 
Eq. (16) with Na = 2. The solutions for these equations 
can be readily obtained. After inserting these solutions 
into Eqs. (18) and (19), we can get the photon spectra 
when t —>■ 00 which are given by 


Pfk — PskiO) 
f^Sk 


—A 6 k‘^Vg/T'^ 


(1 — 2iSkvg/T)‘^ — 6 ^*^“ ’ 


(37) 


Pski, 0 )e^^’^^^ 


(1 -t e^^fe°)(l - 2iSkvg/T) - 2e 

(1 — 2i5kvg/TY — 


2ika 


(38) 


where k = ka + 6 k. 

Before discussing the detailed results, we first look 
at what has been modified by the effective collective 
coupling. Since aj{t) are the one-photon excitation 
amplitudes, we have d{t) = ai{t)\eg) + a 2 {t)\ge) and 
the coupling matrix can be written as V 2 {Sk = 0 ) = 
-(r/‘2){\eg){eg\ + \ge){ge\ + d'^-‘^\eg){ge\+d^<‘<^\ge){eg\). 
It is noted that this coupling matrix is non-Hermitian 
because it includes the effects of dissipation to the pho¬ 
ton modes [i^. The eigenvalues and the eigenstates of 
V 2 {Sk = 0 ) can be calculated to be 

r r 

A± = --(1 ± cos(/caa)) ± i-sin(fcaa), (39) 
|±) = -^{\eg)±\ge)). (40) 

The real parts of the eigenvalues give the collective cou¬ 
pling strength while the imaginary parts give the energy 
shifts. The eigenstates of the coupling matrix are two 
Dicke states. We can change the atomic separation to 
tune the collective couping strength and the energy shift. 
In the following we compare the results for four different 
atomic separations, i.e., a = A/4, A/2, A /8 and 3A/8. 

The results are shown in Fig. 3 where we assume that 
the center of the input pulse is initially 10 /A away from 
the first atom. In Fig. 3(a-c) we present the results when 


a = A/2 and F = Avg. In this case, the eigenvalues of the 
coupling matrix are A+ = 0 and A_ = —F. The collective 
coupling between the two atoms when a = A/2 modifies 
the collective couping strength but does not shift the en¬ 
ergy. The two single-photon excitation eigenstates are 
degenerate. Different from the usual case, here |A_) is 
superradiant while |A+) does not couple to the waveg¬ 
uide modes. This is because of the extra phase factor 
caused by the spatial separation between the two atoms. 
Therefore, the system is prepared in the |A_) eigenmode 
where the two atoms can have almost equal probabilities 
to be excited (Fig. 3(a)). Since only one eigenmode cou¬ 
ples to the ground state, the photon pulse shape (Fig. 
3(b)) and photon spectrum (Fig. 3(c)) before and after 
the interaction are similar to the one-atom cases. How¬ 
ever, the reflectivity here is about 84% which is larger 
than the one-atom case due to the enhanced collective 
coupling strength. 

For a = A/4, the eigenvalues of the coupling matrix are 
A± = —F/2 ± iF/2. The collective interaction only cause 
energy shift but does not modify the collective coupling 
strength. The two eigenmodes have different energies but 
they couple to the ground state with the same coupling 
strength. In Fig. 3(d) we show the excitation probabil¬ 
ities of the two atoms as a function of time where the 
black solid line is the excitation for the first atom while 
the red dotted line is the excitation for the second atom. 
Due to the interference between the two excitation chan¬ 
nels, the first atom has a maximum excitation probability 
larger than 50% which is not allowed for single atom case. 
The second atom has much smaller excitation probability 
and experiences two cycles of excitation and deexcitation 
and reaches a minimum excitation in between due to the 
destructive interference. The photon pulses before and 
after the scattering are shown in Fig. 3(e) where we can 
see that the reflected pulse is similar to incoming pulse 
while the transmitted pulse has three peaks instead of 
two peaks as in the one-atom scattering. The spectrum 
is shown in Fig. 3(f). Different from the single atom 
case, there is a frequency bandwidth around the resonant 
frequency that can not be transmitted. The bandwidth 
where less than 1% transmittance is about 30% of input 
pulse width. This is similar to the photonic bandgap ef¬ 
fects. In principle, we can tune the coupling strength to 
tune this bandwidth. 

When a = A/ 8 , A± = -(F/2)[(I ± v^/2) T iv^/2] 
where the collective coupling affects both the collective 
decay rate and the energy shift. The excitation proba¬ 
bility as a function of time when a = A /8 is shown in 
Fig. 3(g) where we see that the second atom also ex¬ 
periences a population oscillation due to the destructive 
interference. The dip in between does not touch to zero 
because of the unequal excitation of the two eigenmodes. 
The reflected and transmitted spectra is shown in Fig. 
3(h) where the spectrum is very different from the pre¬ 
vious cases. Both the reflected and transmitted spec¬ 
tra are asymmetric. Most of the higher frequency com¬ 
ponents are reflected while most of the lower frequency 
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FIG. 3: (Color online) (a, d, g) The atomic excitations (Al: atom 1 and A2: atom 2) as a function of time. When t = Q the 
center of the input pulse is 10/A away from the first atom. Here, Al is for atom 1 and A2 is for atom 2. (b, e) The pulse 
shapes of the incoming (I, at t = B/T), reflected (R, at t = 20/r) and transmitted (T, at t = 20/r) photon, (c, f, h, i) The 
spectrum (arb. units) of the incoming (I), reflected (R) and transmitted (T) photon. Parameters: T = Aug, (a-c) a = A/2, 
(d-f) a = A/4, (g, h) a = A/8, (i) a = 3A/8. 





FIG. 4: (Color online) (a) The reflectivity as a function of atomic distance with three different coupling strengths: F = O.SAiig 
(blue line with up triangles), F = Avg (green line with squares), and F = 2Avg (red line with circles), (b) The reflectivity 
as a function of frequency detuning with four different atomic separations when F = Avg'. a — Xj2 (black line with asterisk), 
a = XjA (red line with circles), o = A/8 (blue line with up triangles), and a = 3A/8 (green line with down triangles), (c) The 
reflectivity as a function of F /Avg with three different atomic separations: a = A/2 (black line with square), a = A/4 (red line 
with circles), and a = A/8 (blue line with up triangles). 


components can transmit. This can be explained by the 
coupling matrix. The single-photon excitation states are 
split into two states where |-|-) is shifted up while |—) 
is shifted down. Meanwhile, the |-|-) state has a larger 
coupling strength. Therefore the higher frequency com¬ 
ponents have a larger reflectivity. On the contrary, the 
|—) state has smaller coupling strength which result in 


a smaller reflectivity. In addition, there is a dip in the 
reflected spectrum due to the spontaneous emission can¬ 
cellation between the two decay channel and it is a kind of 
Fano-like interference The dip occurs at the po¬ 

sition when Sk = T/2vg. We can also reflect the lower fre¬ 
quency components but transmit the higher frequencies 
by changing the atomic separation such that a = 3A/8. 
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FIG. 5: (Color online) The concurrences as a fnnction of time 
for three different atomic separations (a = 0.5A: black solid 
line; a = 0.25A: red dotted line; a = 0.125A: bine dashed 
line). When t = 0 the center of the input pulse is 10/A away 
from the first atom. Parameter: F = Aug. 


The spectra is shown in Fig. 3(i) where we can see that 
it is similar to the result of a = A/8 but the spectrum 
is opposite. This tunability may find useful as a single 
photon frequency filter. 

The reflectivity as a function of atomic distance is 
shown in Fig. 4(a) for three different coupling strengths. 
We can see that the reflectivity periodically changes with 
the atomic separation. With small coupling strength 
F = 0.5Az;g (blue line with up triangle), the reflectiv¬ 
ity is maximal when the atomic distance is half inte¬ 
ger of the wavelength. However, when we increases the 
coupling strength, the reflectivity increases significantly 
when the atomic separation is one quarter wavelength. 
When F = 2Avg, the reflectivity with a = A/4 is 97% 
which is even larger than that with a = A/2 (95%). We 
also study how the reflectivity changes when the center 
frequency of the Gaussian pulse deviates from the reso¬ 
nant frequency. The results are shown in Fig. 4(b) where 
we consider four atomic separations. When a = A/4 and 
a = A/2, the reflectivity is maximum if the center fre¬ 
quency of the Gaussian pulse is the same as the resonant 
frequency. The full width at half maximum of the re¬ 
flection curve when a — A/2 is about 2.6Avg which is 
larger than that when a = A/4. It is also larger than 
the result in single atom case which is the result of the 
superradiant effect. However, when a = A/8 (a = 3A/8) 
the reflectivity is maximum if center frequency of the 
Gaussian pulse is blue (red) detuned from the resonant 
frequency which is due to the asymmetric coupling of 
the two eigenmodes. The reflectivity as a function of 
coupling strength for three different atomic separations 
are shown in Fig. 4(c). We can see that in all cases the 
reflectivity increases rapidly as we increase the coupling 
strength. The reflectivity can approach 100% when F is 
large comparing with Avg. 

Due to the collective interaction between the two 
atoms, quantum entanglement can also be generated in 
this process. After tracing over the photonic parts in Eq. 
(2), the bipartite quantum entanglement is calculated to 


C{t) = max{0, V2|ai(t)||a2(t)|} (41) 

where C(t) denotes the quantum concurrence of the two- 
atom system at time t [i^. The concurrences as a 
function of time for three different atomic separations 
(a = 0.5A, 0.25A, 0.125A) are shown in Fig. 5. For 
all three cases, quantum entanglement between the two 
atoms can be generated when the photon propagate 
through which is consistent with the results shown in 
[ 4 ^ . When a = A/2, the concurrence increases to a max¬ 
imum value of about 0.17 when the atoms are excited 
and stays at the maximum values for about 2/F before 
decreasing to zero (black solid curve). More interest¬ 
ingly, when a = A/4 the quantum entanglement under¬ 
goes several oscillations (sudden death and revival, red 
dotted line) due to the interference between the two de¬ 
cay channels. When a = A/8, the quantum entanglement 
also undergoes small oscillations and it lasts longer than 
the other two cases (blue dashed line) because one of its 
eigenstates is sub-radiant state. 


V. N-ATOM CHAIN 

The dynamics of the atomic excitations for N-atom 
system with Gaussian input are given by Eq. ( 8 ) with 
j = 1, 2, • • • , Na- The analytical solution for many-atom 
system is difficult to obtain. We need to apply the nu¬ 
merical method to study the dynamics of the system. 
The spectrum of the photon after the interaction can be 
calculated by Eqs. (12)-(19). 

The effective N-atom coupling matrix is given by 

^ikaCL ^ ^ ^ — 

... g^(-^a~2)fca,a 

^i{Na — l)ka.a ^i{Na — 2)kaa _ _ ^ 

For general cases, this matrix can not be analytically 
diagonalized. However, when a = A/2, the matrix can 
be diagonalized and the eigenvalues are given by Ai = 
NaF /2 and X 2 ,...,n^ = 0. Only one eigenmode can couple 
to the waveguide modes and the coupling strength is Na 
times the single atom coupling strength which is the col¬ 
lective superradiant effect. The superradiant eigenstate 
is given by (l/y^) 7 / ^° ( -l)lff, ■.. , 6 ^,... ,. 9 ) which is 
the timed Dicked state [^. For general cases, multiple 
eigenstates can couple to the ground state and they may 
interfere with each other. 

Let us take Na = 5 as an example to study the be¬ 
haviors of the photon transport in this multiple-atom 
system. The results when a = X/2 and F = Avg are 
shown in Fig. 5(a-c) where (a) is the excitation proba¬ 
bility for each atom as a function time when the photon 
passes through, (b) is the photon pulse before and af¬ 
ter the interaction and (c) is the corresponding spectra. 
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FIG. 6: (Color online) (a,d,g) The atomic excitation probabilities as a function of time when the photon passes through. When 
f = 0 the center of the input pulse is 10/A away from the first atom. Here Ai is for the ith atom. (b,e,h) The pulse shapes of 
the incoming (I, at f = O/T), reflected (R, at t = SO/T) and transmitted (T, at t = SO/F) photon. (c,f,i) The spectrum (arb. 
units) of the incoming (I), reflected (R) and transmitted (T) photon. Parameters: F = Aug. (a-c) a = A/2, (d-f) a = A/4, (g-i) 
a = A/8. 




FIG. 7: (Color online) (a) The reflectivity as a function of center frequency deviation for different atomic separations: a = A/2 
(black line with asterisk), a = XjA (red line with circles), a = A/8 (blue line with up triangles), and a = 3A/8 (green line with 
down triangles), (b) The reflectivity as a function of atom numbers for two atomic separations: o = A/2 (red line with squares) 
and a = A/4 (blue line with circles). Parameters: F = Aug. 


In this case, only the superradiant eigenstate is coupled 
with the guided photon and we can prepare the system in 
a timed Dicke state where each atom has an equal prob¬ 
ability to be excited. From Fig. 5(a), we can see that all 
of the atoms are indeed excited and deexcited at almost 
the same rate. The small differences in each amplitude 
come from the time-retarded effect. The reflected and 


transmitted photon pulse shape and spectrum are simi¬ 
lar to the single atom case because the collective coupling 
only affects the collective decay but not the energy shift 
when a = A/2. Since the collective coupling and decay 
rate here is 5 times larger than that of the single atom 
case, the reflectivity here is much larger which is about 
96% comparing with 66% in the single atom case. 
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The results when a = A/4 and T = Avg are shown in 
Fig. 5(d-f). The atomic excitation probability as a func¬ 
tion of time is shown in Fig. 5(d) where we can see that 
the five atoms have quite different excitation dynamics 
and they oscillate with time due to the interferences be¬ 
tween each coupling channel. In addition, the atomic 
excitations last much longer than one atom case because 
some photon energy are trapped in the subradiant states 
of the system. The reflected and transmitted pulses are 
also quite different here (Fig. 5(e)). Both the reflected 
and transmitted pulses have multiple peaks and valleys 
due to the interferences of the reflected and transmit¬ 
ted photon amplitude. A finite bandwidth of frequencies 
around the resonant frequency are almost completely re¬ 
flected which is the demonstration of photonic bandgap 
effect (Fig. 5(e)). The reflection and transmission win¬ 
dow here is about O.SAvg which is larger than that in 
two-atom case and it can also be tuned by changing the 
coupling strength which may be used as a single photon 
frequency filter. Different from the two-atom cases, the 
reflection spectrum here has two dips due to the Fano- 
like interferences between the multiple decay channels 

Elliilil. 


The results when a = A /8 and F = Avg are shown 
in Fig. 5(g-i). Similar to the case when a = A/4, the 
atomic excitation can exist for a very long time and the 
excitation probability can oscillate with time (Fig. 5(g)). 
However, the oscillation period here is longer than that 
when a = A/4. The reflected and transmitted photon 
shapes are shown in Fig. 5(h). The spectra are shown in 
Fig. 5(i) where we see that most of the higher frequency 
components are reflected and most of the lower frequency 
components can transmit. This property may be used as 
a single photon frequency filter. There are multiple dips 
in the reflected spectrum due to Fano-like interferences 
between multiple eigenstates [ 3 , |3113 ■ 

The reflectivity as a function of center frequency devi¬ 
ation for four different atomic separations are shown in 
Fig. 6 (a). Similar to the two-atom cases, the reflectiv¬ 
ity is maximum if the center frequency of input photon 
is resonant with the atomic transition frequency for the 
cases when a = A/4 and a = A/2, while the reflectivity 
is maximum if the center frequency of input photon is 
blue or red detuned from the resonant frequency for the 
cases when a = A /8 and a = 3A /8 due to asymmetric 
coupling. The full width at half maximum of the reflec¬ 
tivity curve when a = A/2 is about 5.2Avg which is much 
broader than the one-atom and two-atom cases which is 
a superradiant effect of this collective many-body system 

M- 


Finally, we study how reflectivity changes as we in¬ 
crease the atom numbers in the chain. The results are 
shown in Fig. 6 (b) where we consider two atomic sepa¬ 
rations, i.e., a = A/2 and a = A/4. When a = A/2, the 
reflectivity increases as the atom number increases and 
the reflectivity can approach 100 % if the atom number 
is large enough. However, when a = A/4 the reflectivity 
does not significantly change when we increase the atom 


number. Therefore, to effectively reflect a photon, one 
way is by increasing the coupling strength and the other 
way is by increasing the collective coupling. To increase 
the collective coupling, we can increase the atom num¬ 
bers but the separation between atoms need to be half 
integer of the wavelength. 


VI. EFFECTS OF THE SPONTANEOUS DECAY 


In previous sections, we have neglected the effects of 
the spontaneous decay to the free space. In this section 
we numerically study how the spontaneous decay of the 
atom to the free space affects the scattering properties. 
Let us take N = 2 atom for example. For two-atom 
system, Eq. (14) including 7 becomes 

- i6kvgx{Sk) = b{Sk) - |x(dfc) -h V 2 {Sk)x{Sk) (43) 

with X = (xiiX 2 )^, b = ibi,b 2 )^. It is readily to solve 
for xi and X 2 from Eq. (43). The spectra after the 
interaction can then be calculate from Eqs. (18) and 
(19) and they are given by 


Psk 


7^ - 4fdfcug7 - ASk'^v^ 

(r -t 7 - 2i5kvgY - r2e2*fc“ ’ 

. (1 + + 7 - 2»dfcug) - 2re 

(F -|- 7 — 2i6kvgY — r 2 e 2 *^“ 


(44) 

2ika 


(45) 


where k = ka + 5k. 

The results are shown in Fig. 8 (a) and (b) where the 
black solid line is the spectrum of the input photon, the 
red dashed line is the spectrum of the reflected photon, 
the blue dotted line is the spectrum of the transmitted 
photon, and the green solid line is the summation of 
the reflected and transmitted photon in the waveguide. 
When there are only two atoms, the result for 7 = r/5 is 
shown in Fig. 8 (a) and the result if 7 = F is shown Fig. 
8 (b). In both cases, the total probability that the pho¬ 
ton stays in the waveguide is reduced because the photon 
can leak out to the free space. However, we notice that 
the general features of the scattering spectrum do not 
change for smaller decay rate when 7 = r/5, i.e., there is 
a range of spectrum that can not transmit. In this case 
the photon has a probability about 68 % staying in the 
guided modes. For larger decay rate when 7 = F, only 
about 26% probability the photon can stay in the waveg¬ 
uide and the general features change where the photonic 
bandgap effect disappears. This is because the sponta¬ 
neous decay can eliminate the interference of the emitted 
photon from the two atoms. However, if we increase the 
atom number such that Na = 5 the photonic bandgap 
effect appears again even if 7 = F (Fig. 8 (c)). This is 
because the collective coupling can be much larger than 
the spontaneous decay rate. 
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FIG. 8 : (Color online) The spectrum (arb. units) of the inpnt (I), reflected (R) and transmitted (T) photon with spontaneous 
decay to the free space. The black solid line is the inpnt spectrum. The red dashed line is the spectrum of the reflected photon. 
The blue dotted line is the spectrum of the transmitted photon. The green dash-dotted line is the total spectrum of the photon 
remains in the waveguide, (a) = 2, T = Aug, and 7 = T/5; (b) = 2, T = Aug, and 7 = F; (c) A'^ = 5, T = Awg, and 7 = F. 


Here, we should note that we have neglected the dipole- 
dipole interaction between the atoms induced by the vac¬ 
uum modes in the free space. This is a good approx¬ 
imation if the atomic separation is not very small (e. 
g. larger than A/10) [^, [^. However, if the atomic 
separation is very small (e. g., smaller than A/10), the 
dipole-dipole interaction between the atoms induced by 
the vacuum modes can be comparable with or even larger 
than the spontaneous decay rate 7 and it should be in¬ 
cluded in the calculations. In this paper, we focus on 
the dipole-dipole interaction between the atoms induced 
by the guided modes which is a long-range effect and the 
smallest atomic separation used in the numerical calcula¬ 
tions is A /8 where the dipole-dipole interaction induced 
by the vacuum modes in the free space can be safely ne¬ 
glected. 


VII. SUMMARY 

In summary, we derive a time-dependent theory to 
study the dynamics of a single photon pulse propagates 
through a linear atomic chain which is coupled to a ID 
single mode photonic waveguide. This theory allows us to 
study the real time evolution of a single photon transport 
and atomic excitations. The shape of the single photon 
pulses and the arrangements of the atoms can be both 
arbitrary. Our result is consistent with the previous cal¬ 
culations when there is only one atom. For more atoms, 
the collective interactions can significantly modify the 
dynamics of the system. The reflectivity can approach 
100 % for a photon with finite bandwidth by changing 
the ratio of coupling strength and the photon linewidth 
or by changing the number of atoms in the chain to in¬ 
crease the collective coupling strength. The spectrum of 
the reflected and transmitted photon can also be signif¬ 
icantly different from that of the single atom case. We 
numerically show the photonic band gap effects where 
a bandwidth of frequency can not transmit and we can 
also tune this gap by changing T /Avg or by changing 
the number of atoms. We can also reflect certain part of 
the frequency components but transmit the other part of 


frequency by controlling the collective coupling between 
the atoms. This may find important applications for sin¬ 
gle photon frequency filter. The quantum entanglement 
between distant atoms can also be created by the guided 
photon which may find important applications in quan¬ 
tum information and quantum computing. Besides, we 
also demonstrate the Fano-type interference and the sup- 
perradiant effects in this system. 

The spontaneous decay to the free space can reduce 
the photon amplitude that stays in the photonic waveg¬ 
uide. If the coupling strength between the atom and the 
waveguide mode is much larger than the spontaneous de¬ 
cay rate, the general feature is similar to the case when 
we neglect the spontaneous decay. However, when the 
spontaneous decay rate is comparable to the coupling 
strength between the atom and the waveguide mode, the 
feature of the reflected and transmitted spectrum can be 
quite different. However, if we increase the atom number 
to increase the collective coupling strength, the reflected 
and transmitted spectrum can have feature similar to the 
case without spontaneous decay rate. 
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APPENDIX: CONCURRENCE OF TWO-QUBIT 
MIXED STATE 

The quantum entanglement of a two-qubit system can 
be well quantified by quantum concurrence which is given 
by m 

C{p) = max{0, Ai — A 2 — A 3 — A 4 } (46) 

where Xi{i = 1,2, 3,4) are the eigenvalues, in decreasing 
order, of the Hermitian Matrix A = y/pPy/p with p = 
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{(Ty(^ay)p*{ay(^ay). 

For the two-atom system in this paper, we have 


The density matrix of the atomic system is pA(t) 
TrpipAFit)) which can be calculated to be 


PAFit) = |'I'W)(^(<)| (47) 


with 

2 

14^(0) = E 0) -t- E U). 

j=l k 

(48) 


PA{t) 


rEJ/3fc(t)p 0 0 

0 |Q;i(t)P ai{t)a2{t) 

0 a{{t)a2{t) \a2{t)\^ 

0 0 0 


O' 

0 

0 

0 . 


• (49) 


The Hermitian matrix A(t) is then given by 


J 


m 


■0 0 

0 \a,ma 2 {tr 

0 \ai{t)\\a2(t)\al{t)a2(t) 
0 0 


0 01 


\ai{t)\\a2it)\ai{t)a2{t) 0 

0 0 . 


and its eigenvalues are given by Ai(t) = which is the Eq. (41). 
^/\ai{t)\\a2{t)\, X2it) = V2\ai{t)\\a2{t)\,X3{t) = 

Xiit) = 0. Hence, the concurrence of the two-atom 
system is given by 


C(t) = max{0, \/|ai(t)||a 2 (i)| - v^|Q;i(t)||a 2 (i)|} (51) 
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